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We present phase space distributions and multiplicities of , K~ and mesons produced in Ar+KCl 
reactions at a kinetic beam energy of 1.756 AGeV and measured with the HADES spectrometer. The inverse 
slope parameters and yields of kaons supplement the systematics of previous measurements. The percentage of 
K~ mesons coming from decay is found to be 18 ± 7%. 

PACS numbers: 25.75.Dw,25.75.-q 



I. INTRODUCTION 



The systematic study of sub-threshold produc- 
tion yields, phase space distributions and flow ob- 
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servables in relativistic heavy-ion collisions at various 
beam energies and for various system sizes and cen- 
tralities has attracted much attention, in particular in 
the context of in-medium properties of A'+ and 
mesons [1-10]. Corresponding experiments have been 
flie focus of studies by the KaoS [11-15] and FOPI 
[17-21] collaborations over the last two decades, fol- 
lowing up on the pioneering experiments at Bevalac 
[22]. The heavy-ion data are supplemented by proton- 
induced reactions on nuclei [23, 24] with measure- 
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ments of inclusive and K production. Con- 
cerning the elementary proton-proton reactions, the re- 
cent measurements by COSY-11 [25], TOF [26] and 
ANKE [27, 28] coUaborations provided valuable ref- 
erence data. 

High statistics data for K'^ production allowed to 
draw conclusions on the /r+-nucleon potential and on 
the nuclear equation of state [11, 13]. From the anal- 
ysis of the out-of-plane [29, 30] and sideward flow 
[20, 31], the ii'+-nucleon potential is deduced to be 
mildly repulsive, and the measured production 
yields are in agreement with the scenario of a soft 
nuclear equation of state (incompressibility modulus 
Kn « 200 MeV) [32, 33]. The interpretation of the 
K~ data is still undergoing systematic evaluations to 
pin down the value of the presumed attractive K~- 
nucleon potential and its momentum dependence. 
Recent data [15] gave an important contribution to- 
wards understanding the production mechanism of 

mesons. Indeed, the combined analysis of 
and K~ suggests that a substantial part of the observed 
K'^ mesons is due to a strangeness exchange mecha- 
nism [34]. Furthermore and K~ exhibit distinc- 
tively different in-medium properties, both in theoreti- 
cal approaches [1-10, 32-35] and in experiments [11- 
13, 15, 19, 20, 23]. 

The spectrometer HADES [36], primarily designed 
to measure di-electrons [37], has been recently em- 
ployed for the identification of strange mesons as well, 
showing high purity and efficiency for particle identi- 
fication and excellent reconstruction capability of sec- 
ondary decay vertices [38]. For the first time, a com- 
bined and inclusive identification of K~^, K~ and (f) 
mesons was carried out in the same experimental set- 
up at a sub-threshold beam energy. ("Sub-threshold" 
refers to free nucleon-nucleon colUsions and applies 
here to the K~ and (/> channels). Our measurement al- 
lows firm conclusions on the fraction of K~ mesons 
that originate from the </> decay. It was argued in [39] 
that, if a substantial part of the observed K~ stems 
from (j) decays, the interpretation of data in terms of 
a strongly reduced effective anti-kaon mass in the nu- 
clear medium needs a proper account of this channel. 
Indeed, the data in [21], unfortunately hampered by 
poor statistics, gave a first hint to the importance of 
K~ feeding from decays. Here, we report on the 
the phase-space distributions and yields of and (j) 
mesons in the reaction Ar + KCl at a kinetic beam en- 
ergy of 1.756 AGeV. We provide an estimate of the 
feeding contribution of </> decays to the K~ abundance. 
Our paper is organized as follows. Section II de- 
scribes the experiment and section IE is devoted to the 
data analysis. Experimental results on transverse mass 
spectra, rapidity distributions and multiplicities of K"^ 
and 4> mesons are presented in section IV. Our findings 
are discussed in section V and summarized in section 
VI. 



n. THE EXPERIMENT 

The experiment was performed with the High 
Acceptance Di-Electron Spectrometer HADES at the 
heavy-ion synchrotron SIS at GSI Helmholtzzentrum 
fiir Schwerionenforschung in Darmstadt, Germany. A 
detailed description of the spectrometer is presented in 
[36]. In the following we summarize the main features 
of the apparatus. HADES consists of a 6-coil toroidal 
magnet centered on the beam axis and six identical de- 
tection sections located between the coils and covering 
polar angles from 18° to 85°. In the measurement pre- 
sented here, the six sectors comprised a gaseous Ring- 
Imaging Cherenkov (RICH) detector, four planes of 
Multi-wire Drift Chambers (MDCs) for track recon- 
struction, in the outermost plane, two opposite MDCs 
were not installed, and they have not been used for 
identification and two Time-of-Flight walls (TOF 
and TOFino), supplemented at forward polar angles 
with Pre-Shower chambers. For each sector, the TOF 
and TOFino/Pre-Shower detectors are combined into 
META (Multiplicity and Electron Trigger Array) de- 
tectors. 

An is_|_Ar beam of ~ 10^ particles/s was incident 
on a four-fold segmented KCl target with a total thick- 
ness corresponding to 3.3% interaction length. A fast 
diamond start detector located upstream of the tar- 
get was used to determine the interaction time. The 
data readout was started by a first-level trigger (LVLl) 
decision, requiring a observed charged-particle mul- 
tiplicity MUL > 18 in the TOF/TOFino detectors, 
accepting approximately 35% of the nuclear reaction 
cross section. It was followed by a second-level trig- 
ger (LVL2) requesting at least one hit in the RICH 
detector per event. This selection has, however, not in- 
troduced any bias on the hadron analysis. A fraction of 
10% of the LVLl triggered events was stored regard- 
less of the LVL2 condition. Comparison between the 
pure LVLl and the LVL1+LVL2 events showed only 
insignificant differences in the total track multiplicity 
and hence the centrality of the reaction. 

III. DATA ANALYSIS 

A. Track reconstruction 

A reconstructed hadron track in the spectrometer is 
composed of inner and outer track segments in the 
MDCs which define straight lines between the first 
two MDC planes, the last two MDC planes behind 
the magnet, and a hit point in one of the META de- 
tectors. In the so-called cluster finder [36, 44] and 
track-segment fitting procedures possible trajectories 
through the MDC segments are calculated. The point- 
ing vector of the outer track segment was used for 
matching the META hit to the MDC track segments, 
generating a track candidate. A Runge-Kutta algo- 
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rithm was implemented for calculating the momentum 
of each of the track candidates by solving the equa- 
tion of motion inside the magnetic field region [45]. 
The quality of the META hit-matching procedure and 
of the Runge-Kutta fitting (characterized by values) 
are used to create a list of ordered track candidates. To 
resolve matching ambiguities, only the track candidate 
with the smallest value of the product of the two ^ 
values is declared a true track. Its components and 
associated track candidates are deleted from the can- 
didate list. This procedure is repeated until no more 
track candidates are left in the list. 



B. Particle identification 

Since the RICH detector is only used to distin- 
guish from hadrons, this system is not employed 
in the following analysis which deals exclusively with 
hadron reconstruction. Hence, only the time-of-flight 
walls TOP and TOFinoH-Pre-Shower (this detector is 
used to refine the position of the TOFino hits) and the 
MDCs are used to identify hadrons. 

The particles can be distinguished by their velocity 
and by their energy loss {dE/dx). The particle veloc- 
ity (3 (in units of velocity of light c) is determined by 
the time-of-flight between the diamond start detector 
and the corresponding hit in one of the two time-of- 
flight detectors and the reconstructed flight path. Fig. 
1 (top panel) shows the velocity j3 of the particles as a 
function of their polarity times momentum for the TOF 
system. Together with the experimental data, curves 
exploiting the relation p = p-fnioc, with mo as parti- 
cle rest mass in vacuum and 7 = (1 — /9^)^^/^, are 
shown for different particle species. The correspond- 
ing distribution for the TOFino system is not shown. It 
has a lower spatial and temporal resolution. 

In addition to the arrival time, also the signal height 
is measured in TOF and TOFino. The latter one is pro- 
portional to the energy loss in the scintillators. The 
so obtained distribution of the particle energy loss in 
the scintillator of the TOF system is shown in Fig. 1 
(TOF dE I dx, center panel) as a function of polarity 
times momentum. The solid lines are obtained by fit- 
ting the Bethe-Bloch formula [42] to the proton distri- 
bution and then scaling the other curves according to 
the different particle masses. The separation power of 
this observable is comparable to the one obtained from 
the velocity measurement. 

The bottom panel of Fig. 1 shows a similar dis- 
tribution, but for the combined energy loss in the 4 
MDCs of one sector. The energy-loss information is 
derived from the sense wire signals after unfolding ef- 
fects due to the arrival time distribution of the drifting 
electrons which vary with track incident angle [38, 44]. 
The solid curves in the same panel were calculated us- 
ing the Bethe-Bloch formula. For the identification 
of TT^ and protons detected in the TOF, the /3-versus- 
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FIG. 1: Color online. Particle velocity /3 (top, for the TOF 
detector) and dE/dx (middle for TOF detector, bottom for 
the MDC) as a function of the momentum times the particle's 
polarity sign for the TOF detector. 



momentum information is used and the particles are 
selected by applying a graphical cut around the solid 
lines shown in Fig. 1 (top). The identification of tt^ 
and protons in the TOFino is performed in the same 
way, if the hit multiplicity per paddle is equal to 1; 
in case of multiple hits, only the coiTesponding MDC 
dE/dx information is used to perform the selection. 

For the identification, only the TOF data are 
considered, because of rather large background in the 
TOFino due to its limited granularity and time reso- 
lution. As evident from Figs. 1 (top) and 2, the huge 
vr^ and proton background hides completely the 
signals. In order to reduce this background, graphi- 
cal cuts based on calculated dE/dx distributions have 
been applied. The effect of these cuts on the signal- 
to-background ratio for the signal is exhibited in 
Fig. 2, where the application of the TOF dE/dx and 
the MDC dE/dx vs. momentum selection reduces the 
background such that and signals are clearly 
visible above the background. 



C. K+ and K spectra 

Examples of the mass distribution obtained by ap- 
plying both the TOF dE/dx and MDC dE/dx cuts 
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FIG. 2: Color online. Polarity times mass over charge distri- 
butions for particle tracks measured in the TOF detector for 
different cut conditions. The spectrum on top (blue curve) 
shows the original distribution, whereas the lower two spec- 
tra show the effects of the graphical cuts in the TOF dE/dx 
(green histogram) and additionally the MDC dE/dx (red 
histogram) distributions, respectively. 
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are shown in Fig. 3. The data are shown in intervals of 
rapidity, y = \ \\i{{E +pz) / [E — p^)), and subtracted 

transverse mass {mt-rnK± with nit — \J + p^), 
where pt, Pz and E are the transverse and longitudi- 
nal momenta and the total energy of the particles, re- 
spectively. The solid lines in Fig. 3 corresponds to 
a combined fit of the signal and background, while 
the dashed lines indicate the background. A Gaus- 
sian fit of the signal after the background subtraction 
and integrating over the whole covered phase space 
gives the following mean values and dispersions for 
the reconstructed masses: {niii±) — 485MeV/c^, 
(Tif± = 20 MeV/c . The slight deviation of about 2% 
from the nominal rest mass of = 494MeV/c^ 
[42] is attributed to imperfections in the time-of- 
flight calibration, to be improved for further analysis. 
Note, however, that all kinetic quantities are derived 
from the measured momentum while the time-of-flight 
serves for particle identification only. The signal-to- 
background ratio (S/B) for kaons depends on the lo- 
cation in phase space. It varies for mesons from 
1.1-35.6 and for mesons from 0.5-4.1. 

The raw phase space distribution for the identified 
and mesons after background subtraction is 
exhibited in Fig. 4, where the color code refers to 
the counting rate. Due to low statistics or large cor- 
rection factors, mainly connected to the detector ac- 
ceptance, not all of the shown bins are used in the 
further data analysis: only the rapidity bins between 
0.1 < ulab < 0.7 have been considered. The three 
dashed lines in Fig. 4 refer to the laboratory polar an- 
gles of 18°, 45° and 85°, respectively. The accep- 
tance is Umited to the TOF region 45° < 6* < 85°. 



FIG. 3: Reconstructed mass distributions of K (top) and 
mesons (bottom) for different phase space regions (left: 
0.5 < y < 0.6 and WOMeV/c^ < rrit < 150MeV/c^ 
right: 0.5 < y < 0.6 and lOOMeV/c^ < rut < 
150i\f eV/c^). The solid lines represent combined fits to the 
signal and background, whereas the dashed lines show the 
background part only. 



D. Efficiencies 

In order to extract quantitative information on the 
yields and momentum distributions over the whole 
phase space, acceptance and efficiency corrections 
must be applied. The total yields are finally ob- 
tained by extrapolating the measured phase space to 
47r. The geometrical acceptance for and 
has been determined as a function of the emission ra- 
pidity and transverse mass. A full-scale simulation 
of the HADES spectrometer was performed with the 
GEANT3 package [46]. Kaons and antikaons were 
generated with a flat (white) distribution in rapidity 
(0 < yLAB < 0.7) and transverse mass (0 MeV/c^ < 
nit — f^K < 450 MeV/c^). Then the generated par- 
ticles were propagated through the detectors and the 
resulting hits subjected to the selection criteria applied 
to the real data. The reconstruction efficiency has been 
evaluated in two steps as a function of y and mt — rriK- 
First, simulated and tracks were embedded in 
real Ar + KCl experimental events and the single-track 
reconstruction efficiency was estimated. In a second 
step the PID efficiencies for were calculated us- 
ing experimental data only, in order to reduce the sys- 
tematic uncertainty introduced by the digitization. The 
efficiency of the TOF dE/dx and MDC dE/dx cuts 
were calculated selecting a kaon sample using only one 
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rapidity y 

FIG. 4: Color online. Distribution of measured K'^ (top) 
and (bottom) yields in the TOF detector as a function of 
the subtracted transverse mass mt — m^ ± and the laboratory 
rapidity y of . The color code refers to the counting rate. The 
dashed lines indicate the polar angles Q = 15°, 45° and 85° 
in the laboratory system. The dotted curves are for constant 
kaon momentum starting from p — 200 MeV/c in steps of 
200 MeV/c. The vertical solid line depicts mid-rapidity. 
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FIG. 5: Color online. Reconstruction efficiency for K and 
A'" as a function of the transverse mass for different rapidity 
intervals. The acceptances are not included (see text). 
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of the two cuts and estimating the reduction of the sig- 
nal in the sample when the second cut is applied [43]. 
This procedure was applied for different rapidity and 
transverse mass intervals. The total reconstruction ef- 
ficiency was obtained by multiplying the three com- 
ponents: track reconstruction, TOF dE/dx cuts and 
MDC dE/dx cuts. Assuming that the efficiency for 
a given momentum slice is independent of the emis- 
sion angle, the total efficiency was extrapolated also 
for the y — nit bins with low statistics. The recon- 
struction efficiency for and K~ has also been es- 
timated using only simulations, exploiting the embed- 
ded tracks method. The digitized data have been tuned 
such that the simulated dE/dx distributions of both 
MDC and TOF detectors reproduce the experimental 
distributions. The so obtained efficiencies show a very 
good agreement (within 5%) with the efficiencies ex- 
tracted from the experimental data and provide the ba- 
sis for a reaUstic estimation of the systematic errors. 
The resulting efficiency is shown in Fig. 5 as a func- 
tion of nit — rriK for different rapidity bins. These 
values do not include the acceptance which varies from 
20-40% for both particles analyzed only in the TOF re- 
gion and in the four sectors equipped with four MDC 
planes each. The phase space coverage is displayed in 
Fig. 4. 



In order to reconstruct </> mesons, K'^ and K~ , iden- 
tified in both the TOF and TOFino system and in all six 
sectors, are combined to pairs after the application of 
quality cuts which are more selective in the case of the 
TOFino detector. 

The resulting invariant-mass distribution for the — 
pairs is shown in Fig. 6 (top) where the (f) signal is 
clearly visible. Despite of the fact that the kaon iden- 
tification is more selective for the particles which hit 
the TOF than the TOFino detector, the presence of the 
peak at the nominal (f) mass with signal to background 
ratio better than 1 (see Fig. 6) shows that the purity 
of our selection is satisfactory for this investigation. 
Nevertheless, the continuum visible under the peak 
is contaminated with misidentified kaons and can- 
not be considered for quantitative conclusions on the 
K^—K^ non-resonant pair yield. To extract the (/> sig- 
nal, the mixed-event technique was employed for the 
determination of the combinatorial background. Since 
a four-fold KCl target stack was used, only events were 
combined for which the reaction took place in the same 
target segment. The selection was done by calculating 
the minimal distance of the global event vertex w.rt. 
the nominal target positions. Additionally, only events 
belonging to the same centrality class {MUL ± 4) 
were combined. The resulting combinatorial back- 
ground distribution is shown by the hatched histogram 
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FIG. 6: Invariant-mass distribution of K'^ — K~ pairs (top). 
Tiie combinatorial background (shaded area) is obtained by 
the mixed-event technique. The baclcground-subtracted dis- 
tribution (bottom) shows a </> meson signal (grey area with a 
Gaussian fit) with a total yield of 168 ± 18 counts. 
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FIG. 7: Color online. Invariant transverse-mass spectra of re- 
constructed A"^ (left) and K~ (right) mesons. For the sake 
of clarity the spectra are plotted for several laboratory rapid- 
ity bins, see legends. The spectra are scaled by the factors as 
indicated in the plots. 



in Fig. 6 (top), where the normalization has been ob- 
tained by scaling the mixed-event distribution to the 
same-event distribution (arithmetic sum) in the invari- 
ant mass region 1050 - 1400 MeV/c^. 

Fig. 6 (bottom) shows the signal distribution after 
the background subtraction. We obtain for the recon- 
structed (j) meson mass = (1017.8 ± 0.9) MeV/c^, 
a width of = (6.2 ± 0.8) MeV/c^ and a total statis- 
tics of 168 ± 18 counts. In order to extract the geo- 
metrical acceptance and the reconstruction efficiency a 
GEANT3 simulation of a "white" meson spectrum 
has been propagated through all the analysis steps. 
Each simulated — pair from a decay has 
been embedded in a real event to provide a realistic 
environment. The inclusive reconstruction efficiency 
for the (f) meson has been evaluated to be about 1%. 
In addition the fraction of 90% of all is lost due to 
the acceptance and furthermore only 49.2% of all 
mesons decay into the channel or , when us- 
ing the free decay branching ratio. 



IV. EXPERIMENTAL RESULTS 

A. Transverse mass spectra 

The invariant transverse mass spectra for , K~ 
and (f) for different laboratory rapidity bins are exhib- 
ited in Figs. 7 and 8, respectively. These distributions 
were obtained by applying the efficiency correction to 
the raw data, including also the decay probability of 
each particle in channels that are not reconstructed in 




FIG. 8: Invariant transverse-mass spectrum of reconstructed 
(j) mesons for the rapidity range 0.2 < y < 0.6. 



this analysis. The distributions show the number of 
counts per LVLl trigger, per transverse mass and ra- 
pidity unit, divided by TOj . This representation is cho- 
sen to ease the comparison with a Boltzmann distribu- 
tion. A fit according to 

~1 — ^=C'(y)exp — — (1) 

dnitdy \ TB{y) J 

yields the solid straight lines shown in Figs. 7 and 
8, which describe indeed the data quite well. The in- 
verse slopes of each distribution, Tsiy), is shown in 
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FIG. 9: Rapidity dependence of the inverse slope parameters 
of K+ (top) and K' mesons (bottom). The full symbols 
show the measured data whereas the open ones are reflected 
data with respect to the cm. rapidity. The lines represents 
fits with Eq. (2). 



Fig. 9 as a function of the center-of-mass (cm.) ra- 
pidity (yc.m. = y — y{c.m.), where y{c.m.) = 0.858 
for symmetric colhsions at 1 .756 AGeV) for (top) 
and K~ mesons (bottom). The full symbols display 
the measured data, whereas the open ones are reflected 
data with respect to the cm. rapidity. The curves in 
Fig. 9 represent fits with 



Toff 
cosh(j/) ' 



(2) 



assuming a thermal source. The obtained parameter 
Teff represents the inverse slope at mid rapidity and 
may be considered as an effective temperature of the 
kinetic freeze-out of the respective particle. In partic- 
ular, one can see from Fig. 9 that Toff is lower for K~ 
than for A'+. This can be attributed to different freeze- 
out conditions for the two meson species [12]. 



B. Rapidity distributions 

The fitted and invariant transverse-mass 
distributions <PN/{dm.tdy) in Fig. 7 are integrated 
within the interval < — rriK < cxo in order 
to obtain the meson yield per rapidity unit. A differ- 
ent approach to determine this integral is to integrate 



the measured data points and use the fits only for the 
extrapolation to the unmeasured phase-space regions. 
However, the difference in yield to the used method is 
negligible. The resulting rapidity density distributions 
of and are displayed in Fig. 10. The full cir- 
cles show the values calculated by the integration of 
Eq. (1), i.e. by means of 



dN 
dy 



2moc^Tl{y,) + 2T%{y, 



(3) 



The parameters Tsiyi) and C{yi) are taken from the 
fits with Eq. (1) to the invariant transverse mass spec- 
tra. The open circles represent the reflection with re- 
spect to the center-of-mass rapidity. The rapidity den- 
sity distributions have been fitted with Gaussian func- 
tions (solid curves represented in the two panels of 
Fig. 10). The distribution for the has a larger width 
than for the K^, a fact which is also observed in other 
measurements [15]. 

In the fitting procedure, the center of the Gaussian was 
fixed at zero, taking into account the symmetry of the 
reaction. By integrating the fit function one obtains 
the kaon multiplicity per triggered event (LVLl). If 
the reflected points are included the extrapolation of 
the measured yield into the unmeasured rapidity region 
amounts to 35% and 31% of the determined total yield 
for if + and K^, respectively. The resulting multiplic- 
ity values are collected in Tab. I. The first error follow- 
ing the multiplicity value indicates the statistical error 
The evaluation of the systematic eiTor is done in two 
steps. First, the graphical cuts on the MDC and TOE 
dE/dx distributions have been varied in width within 
30% and the corresponding variation of the multiplic- 
ity calculated. This results in the first systematic error 
following the statistical one. 

As shown in Fig. 10, the Gaussian parametrization of 
the data reproduces in a satisfactory way the exper- 
imental distribution. To check the sensitivity of the 
multiplicity results to the assumed parametrization, a 
double Gauss fit was also applied to the data following 
the procedure described in [47]. The resulting uncer- 
tainty on the kaon multiplicity is shown by the third 
error in Tab. I. The systematic eiTor of the extracted 
inverse slope values has been estimated by varying the 
graphical cuts on the MDC and TOE dE/dx distribu- 
tion of 30% as already mentioned above. 
Figure 10 shows in addition to the experimental data 
and the Gauss fits the rapidity density distributions ob- 
tained from simulations of an isotropic thermal distri- 
bution of and K~ respectively. The simulated dis- 
tributions are represented by the dotted lines in upper 
and lower panels of Fig. 10 and is normalized to the 
total integral of the experimental ones. The measured 
distributions are systematically wider than the simu- 
lated ones (see Tab. I). For the increase in width 
is larger and amounts to roughly 30%. 
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FIG. 10: Rapidity density distributions of (top) and K~ 
mesons (bottom). The full symbols show the measured data 
whereas the open ones are reflected data with respect to the 
cm. rapidity. The curves represent Gaussian fits to the data. 
The overlayed dotted lines are simulated isotropic thermal 
distributions which have been normalized to the data (see 
text). 



Due to the extended rapidity bin chosen for the <f> 
mesons the change of the yield within the bin is not 
negligible. Instead of using the center rapidity of this 
bin we take a value weighted with the expected thermal 
distribution across the bin width. This procedure re- 
duces the extrapolated total yield by around 24% com- 
pared to an extrapolation where the rapidity point is in 
the center of the bin. The systematical error estimated 
for the total (p multiplicity/LVL 1 shown in Tab. I is ob- 
tained by varying the temperature Toff within its statis- 
tical eiTor. 

Existing theoretical calculations assume an isotropic 
emission of the <j) meson [39-41], but since the angular 
distribution of the emitted cf) has not yet been measured 
at these energies, an additional source of systematical 
error must be considered. On the base of the before 
mentioned maximal deviation of the rapidity dis- 
tribution from an isotropic source, we assume 30% to 
be the maximum error in the evaluation of the width of 
the rapidity distribution for the (j> meson. This trans- 
lates into a second systematic error for the cj) multi- 
plicity/LVLl which is also presented in Tab. I. This 
has been calculated assuming a larger width of the (f> 
rapidity density function and hence leads to an asym- 
metric error. 

The value of the effective temperature of the (p meson 
is shown as well in Tab. I together with its statistical 
error. 



V. DISCUSSION 



Deviations from a thermal picture of a particle emitting 
source are due to anisotropics of the angular and/or of 
the momentum distributions. The limited acceptance 
of the HADES detector precludes such a differentia- 
tion at the current stage of analysis. 

The rapidity density distribution of the measured </> 
mesons contains only one data point. Therefore, one 
has to fix the width of the distribution in order to esti- 
mate the <f> multiplicity. If a thermal description of the 
rapidity density distribution is assumed, its width ay 
can be calculated according to 

= (4) 

where represents the effective inverse slope for (j> 
mesons. This parameter is calculated taking the in- 
verse slope Tb obtained in fitting the invariant trans- 
verse mass distribution of the <f> meson for the rapid- 
ity bin shown in Fig. 8 and substituting this value into 
Eq. (2). 

Table I shows the total (p meson multiplicity per 
triggered event (LVLl) obtained by integrating the so- 
obtained Gaussian curve. The first error shown is again 
purely statistical. 



Let us first compare our data with results recently 
published by the KaoS collaboration [15]. If the 
/K^ multiplicity ratio is plotted as a function of 
the beam energy for the inclusive data and for vari- 
ous systems measured by KaoS [15, 16], as shown in 
Fig. 11, one can see that the ratio increases with beam 
energy and that the HADES data point fits nicely into 
the energy systematics. 

In Fig. 12 the centrality dependence of the /K^ 
ratio is exhibited for the AuH-Au and Ni-nNi systems 
and for two different beam energies. As a measure 
of the centrality, the number of participants Apart is 
chosen. This quantity is calculated within a geomet- 
rical model of interpenetrating spheres where the im- 
pact parameter range is derived from a comparison of 
the charged particle multiplicities, as measured in the 
TOF/TOFino detectors, to that provided by GEANT 
predictions with the UrQMD transport model [49] as 
event generator The full triangle represents the data 
point measured for ArH-KCl at 1.756AGeV in the 
present experiment. For the LVLl trigger we attribute 
to our data a value of Apart — 38.5 ± 4, being about 
twice the minimum bias value. Our data point fits well 
into the systematics when taking into account the beam 
energy dependence. Looking at the empty squares 
(NiH-Ni at 1.5 AGeV) and full squares (AuH-Au at 1.5 
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TABLE I: Total multiplicity/LVLl and effective inverse slope parameters T^e for K'^ , K~ and 0. The first error refers 
to statistical uncertainties, while the second and third errors represent the systematic ones, as described in the text. The 
fitted Gaussian width of the rapidity density distributions (a^xp) as well as the calculated width for a thermal and isotropic 
distribution {(Jtherm) of the corresponding effective temperature T^s are also listed with systematic and statistical errors. 
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FIG. 11: K~ /K'^ ratios as a function of the kinetic beam 
energy for various systems. The present HADES data point 
is depicted by the full triangle. KaoS data are from [15]. 



AGeV), it can be seen that the ratio stays fairly con- 
stant with Apart, supporting our implicit assumption 
on the centrality independence in using an averaged 
value of Apart- Together with the data, / ra- 
tios predicted by the thermal model [52, 53] are shown 
for kinetic beam energies of 1.5, 1.76 and 1.93 AGeV. 
These calculated values are found to be in fair agree- 
ment with the experimental values. The fact that the 
ratio K~ /K^ does not depend noticeably on the cen- 
trality of the reaction and that the meson yield is con- 
sistent with the prediction of thermal models has been 
used in the past to support the hypothesis that most of 
the K~ are produced in multistep processes like 



NN YK+N, 
Yn ^ NK^, 



Eth = 1.58GeY, (5) 
(6) 



where Y stands for the hyperons A and S, and 
Eth denotes the beam energy threshold in free 
nucleon-nucleon collisions. The reaction (6) is called 
strangeness exchange. The dependence of the if + 
and multiplicities on Apart has been found to be 
quite smooth [15]: Mx+ oc Apart with a^+iNi) = 
1.26 ± 0.06 and M^- oc A^art with aK-{Ni) = 
1.25 ± 0.12. (In contrast, the tt^ multiplicity is found 
to linearly dependent on Apart, en = 1). Along these 
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FIG. 12: Apart dependence of the K~ /K^ ratio for differ- 
ent beam energies and colliding systems. The full triangle 
depicts the present HADES data. The other symbols display 
KaoS results from [15]. The dashed lines show the predic- 
tion of a thermal model [52, 53] for kinetic beam energies of 
1.5, 1.76 and 1.93 AGeV (from bottom to top). 



considerations it has been claimed in [52] that reaction 
(6) is the dominant channel for K~ production. 

Further possible channels for the production of K~ 
in nucleon-nucleon collisions are 

NN (t>NN 

K+K-NN, Eth ^ 2.59 GeY, (7) 
NN K+K-NN, Eth = 2A9GeY, (8) 

with thresholds being higher than for the previous re- 
action. 

The threshold for production induced by a A reso- 
nance A^A — *■ K^K^NN depends on the mass (en- 
ergy) of the resonance, but it is comparable to the 
value quoted in Eqs. (7) and (8). Further secondary 
collisions have been considered in theory [39^1] to 
contribute to the 4> production such as ttN — > (f)N, 
7rA^(1520) (t)N, pN (jiN, pA <pN, 7rp -> 
It has to be pointed out [39-41] that the predicted (f) 
multiplicity in sub-threshold heavy-ion collisions, cal- 
culated including all the channels mentioned above, 
stay below the number extracted from the experimental 
data [21]. 

The newly measured HADES data provide for the 
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first time at SIS energies a consistent measurement 
of K'^ , K" and (f> in the same data sample. The 
4)/K~ ratio measured with the HADES spectrometer 
for Ar+KCl at 1.756 AGeV is found to be 0.37 ±0.13, 
see Tab. 1. This value translates into a fraction of 
18 ± 7% of coming from (f) decay. This implies 
that also below threshold the (/> production contributes 
significantly to the rate, as already argued in [39]. 
Taking for the meson the decay path cr w 46 fm 
as follows from the vacuum width, one can estimate 
a fraction of about 80% decaying outside the collision 
zone, for which we employ a radius of about 10 fm. 
The question of the in-medium width of the me- 
son is certainly far from being settled, since the few 
available data do not yet deliver a consistent picture 
[54]. On the other hand, if we assume that a large 
amount of K^K^ decays would take place in 
the fireball volume, due to the strong interaction expe- 
rienced by kaons in medium the resulting signal of in 
the K~K^ invariant-mass spectrum would be washed 
out. Hence, the K~ stemming from (f) decays recon- 
structed in this analysis are produced mostly outside 
the medium. This effect will also dilute any observed 
medium effects of the . In order to extract precise 
information on the equation of state [8-10] and on the 
ii'^-nucleon potential, one should improve the statis- 
tics of such independent measurements for and 0, 
i.e. tag the according to their production mecha- 
nism. This kind of systematic measurements should be 
done for different systems and different beam energies. 

Also the non-resonant A'+ — production (say 
by the reaction Eq. (8)) plays in pH-p reactions an im- 
portant role, contributing to about 50% of the overall 
— yield [28]. This contribution is not yet 
known in heavy-ion collisions, since the K'^ — K~ in- 
variant mass spectrum measured by HADES still con- 
tains a non-negligible contribution of fake candidates 
(see discussion of Fig. (8)). If one assumes that the 
contribution of the non resonant if + — produc- 
tion in ArH-KCl at 1 .756 Ge V is comparable with the 
value measured for pH-p at 2.7 GeV [28], taking the ra- 
tio (p/K^ measured by HADES, it can be estimated 
that about 38% of all produced K~ do not stem from 
strangeness exchange reactions but from the sum of 
resonant and non-resonant — production (Eqs. 
7,8). 

This consideration further deprives the strangeness 
exchange mechanism of its leading role in the pro- 
duction in sub-threshold heavy-ion reactions. 

In Fig. 13, the (t>/K~ ratio measured by HADES is 
compared to the ratios measured at higher energies in 
PbH-Pb [56] and AuH-Au colHsions [55, 57] and in pH-p 
reactions at 2.7 GeV [28]. For heavy ion collisions 
one can see that, while the ratio is rather constant at 
high energies, it is substantially larger in case of the 
HADES measurement. 

The various curves in Fig. 13 show the prediction 
by a statistical model for three different values of the 
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energies [55-57]. The empty circle presents the value mea- 
sured in p-l-p at 2.7 GeV [28]. The lines are predictions of a 
statistical model [53] for three parameters of the correlation 
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correlation length Rc [53]. This parameter is the ra- 
dius of the volume inside which strangeness is locally 
conserved in the calculation. One can see that a value 
of Rc < 3.2 fm is needed to fit the prediction by 
the statistical model to the HADES point. A smaller 
value of Rc translates into a reduced volume for lo- 
cal strangeness production and conservation. The fact 
that the HADES result is consistent with a smaller vol- 
ume could be interpreted as a hint that the contribution 
to production by strangeness exchange, which is 
more probable for larger volumes, is diminished in the 
SIS energy regime. 

The 4)/K~ ratio measured in elementary reactions is 
thought to be consistent with Rc ~ 1 fm. The fact 
that the ratio measured by HADES is smaller than the 
one extracted for pH-p reactions tight above threshold 
shows the role played by secondary processes in the 
production in sub-threshold heavy-ion collisions. 
On the other hand, if the strangeness exchange process 
would be the dominant mechanism for production, 
the (t>/K^ ratio would be much smaller than the mea- 
sured value. 



VI. SUMMARY 

In summary, we report on a first measurement of 
charged kaons and the (f) meson production in the same 
experiment at SIS energies. The HADES data on 
fit well into the systematics of the results published 
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by the KaoS collaboration. The K~ /K'^ ratio is con- 
sistent with the value predicted by the thermal model. 
The effective inverse-slope parameters of and K~ 
are determined; the inverse slope of mesons is 
found to be higher than that for JsT", in agreement with 
the KaoS data. 

The (t)/K~ ratio has been deduced from our data 
with improved accuracy. The HADES value is found 
to be 0.37±0.13, translating into a fraction of 18±7% 
of K~ stemming from </> decays. Since the de- 
cays that can be measured via the invariant-mass re- 
construction of — K~ pairs are essentially those 
happening outside the nuclear medium, this value may 
be considered as a lower Umit. On the other hand, also 
the non resonant —K~ production may play an im- 
portant role for the total K~ yield, but in order to draw 
quantitative conclusions the purity of the kaon identi- 
fication in heavy-ion reactions must be improved. In 
general it can be concluded that for heavy-ion colli- 
sions at beam energies below the free production 
threshold, the K~ production cannot be explained ex- 
clusively by the strangeness exchange mechanism. 
Furthermore, the 4'/K~ ratio for the Ar-nKCl system 
has been compared to the ratio measured in heavy-ion 
collisions at higher energies and to the ratio extracted 
from pH-p reactions at 2.7 GeV. The 4>/K~ ratio mea- 
sured by HADES has been found to be higher than the 
almost constant value found at higher energies. A cor- 
relation parameter Rc < 3.2 fm has to be used in the 
calculation by the Statistical Model to reproduce the 
HADES measurement. This value is lower than the 
one used to fit the (/>/ K~ measured in heavy-ion colli- 
sions at higher energies. A reduced value of Rc can be 
interpreted as a reduced volume in which strangeness 
production and conservation is taking place and hence 
to a less significant contribution of the strangeness ex- 
change mechanism, to the K~ production, which is 
more probable for larger volumes. 
The contribution of the different secondary processes 
to the total K~ is visible in the difference between the 
(p/K^ ratio measured in pH-p reactions at 2.7 GeV and 
in ArH-KCl collisions at 1 .756 AGeV. The observed rel- 
ative high ratio measured by HADES seems to deprive 
the strangeness exchange mechanism from its domi- 
nant contribution io K~ production, as stressed above. 

It is clear that, in order to extract an excitation func- 
tion, further systematic measurements of the and 
production have to be performed. HADES is a suitable 
apparatus for this kind of studies, especially due its ca- 
pability to reconstruct the (j) not only via the charged 
kaon decay but also via the direct electromagnetic de- 
cay into lepton pairs. Due to the fact that leptons do 
not undergo strong interactions and therefore are able 
to leave the reaction zone nearly undisturbed, one will 
be able to distinguish between those mesons which 
decay inside the medium and those who decay outside. 
This will provide an observable directly sensitive to the 
predicted medium modifications. 
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